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ANALYSIS  OF  SLR  TARGETS  FOR  JASON 


1  INTRODUCTION 

This  report  presents  an  analysis  of  a  design  for  a  retroreflector  array  for  the  JASO.X  spacecraft 
and  mission.  The  analysis  includes  the  direct  numerical  computation  of  a  single  retrorefiector's  and 
array’s  optical  properties  as  well  as  the  space  qualification  characterization  of  the  array's  mechanical 
properties.  The  purpose  of  the  array  is  to  provide  maximized  cross  section  return  for  the  Field 
Transportable  Laser  Radar  Station  (FTLRS)  telescope  on  the  island  of  Capraia  Island. 

Satellite  Laser  Ranging  (SLR)  data  is  an  extremely  powerful  tool  for  precise  orbit  determination. 
Direct  detection  SLR  provides  position  estimation  comparable  with  differential  GPS  and  is  used  as 
the  referenced  "truth"  by  the  scientific  community  for  geoscience  and  navigation.  Although  SLR 
data  serves  as  a  high  precision  orbit  estimator,  it  is  weather  dependent.  Therefore,  this  data  type 
is  uniquely  suited  for  independent  system  performance  validation  of  onboard  spacecraft  systems 
and  periodic  calibration. 

The  report  opens  with  a  review  of  SLR  systems  establishing  both  general  terminology  and 
specific  assumptions  relevant  to  JASON  in  Section  2.  Section  3  considers  single  retrorefiectors  and 
Section  4  discusses  the  analysis  of  the  retroreflector  array.  Core  results  from  the  space  qualification 
testing  of  the  retroreflector  array  are  summarized  in  Section  .5. 
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Figure  1:  Schematic  diagram  of  Satellite  Laser  Ranging. 


2  SATELLITE  LASER  RANGING  SYSTEM  CHARACTERISTICS 

Satellite  Laser  Ranging  in  the  context  of  this  report  is  direct-detection  radar  in  the  optical  wave¬ 
length  regime.  When  an  orbit  is  properly  sampled  and  the  SLR  data  acquisition  system  is  properly 
configured  and  calibrated,  residuals  and  accuracies  on  the  order  of  centimeters  are  routinely  ob¬ 
tained  using  this  method. 

Figure  1  illustrates  the  basic  aspects  of  the  technique.  Time-tagged  time-of-flight  differences  are 
recorded,  corrected  for  system  delays,  and  converted  to  ranges.  These  ranges  then  provide  input 
to  an  orbit  determination  model  which  is  used  to  generate  the  three  dimensional  estimate  of  the 
spacecraft’s  orbit  and  position. 


2.1  Return  Pulse  Detection  and  System  Trade-oflfs 

As  discussed  in  Degnan’s  review  [1],  the  number  of  photoelectrons  generated  by  the  return  pulse 
detection  is  given  by  the  laser  radar  link  equation 


The  factors  in  Equation  1  are: 


JASO^■SLR  Taryet  Analysis  vJ.O  XRL  Mejjio  7971 


3 


VD 

Detector  Quantum  Efficiency 

O^LRCS 

Laser  Radar  Cross  Section 

Eo 

Transmit  Energy 

R 

Slant  Range 

X 

Wavelength 

■4/? 

Receiver  Telescope  Area 

h 

Planck's  Constant 

VR 

Receiver  Efficiency 

c 

Speed  of  Light 

Ta 

One-way  Atmospheric  Transmission 

IT 

Gt 

Transmission  Efficiency 
Transmitter  Gain 

Tc 

One-way  Cirrus  Cloud  Transmission 

In  cases  where  the  system  figure  of  merit  is  the  number  of  return  pulse  photons  reaching  the  front 
aperture  of  the  receive  telescope,  A-.teie:  then 


A  7.  tele  — 


A  pe 

mm 

E„  (A)  ,„GV<tlrcs  (j^) 


(2) 


Equation  2  is  an  absolute  number  of  photons  and  depends  on  the  telescope  aperture  Ap.  If  the 
photon  flux  itself  is  needed  for  ground  station  trade  studies,  it  is  given  by 


ApTlpljD 

Since  the  JASON  Specification  in  Section  2.2  use  a  fixed  ground  station,  the  results  of  this  report 
are  given  in  terms  of  A’pe  and  AA  for  a  13  cm  aperture  only. 

The  factors  in  Equation  1  determine  the  SLR  trade-space  and  can  be  grouped  four  categories: 
transfer  efficiencies,  transmitted  pulse  magnitude,  environmental/orbit  parameters,  and  geometric 
factors.  The  three  transfer  efficiencies,  r/£),  rjx,  and  rjp,  are  fixed  by  the  technology  of  the  ground 
station.  The  initial  number  of  photons,  Eq  is  determined  by  the  laser  source.  Three  of  the 

remaining  factors  in  Equation  1  are  environmental  parameters  determined  by  the  satellite's  orbit 
and  the  specific  ground  station  location  and  time  of  the  measurement  {i.e.  R,  Ta,  or  Tc).  These 
environmental/orbital  parameters  are  beyond  the  experimenters’  control.  Two  of  the  geometric 
factors,  6’t  and  Ap,  are  associated  with  transmitting  and  receiving  the  laser  ranging  pulse  at  the 
ground  station.  Assuming  Gaussian  beam  profile,  the  transmitter  gain  is 


(4) 


where  Or  is  the  divergence  half-angle  and  ^point  is  the  pointing  uncertainty.  The  remaining  geometric 
factor,  aLRcs,  is  determined  by  the  SLR.  target  and  will  be  a  primary  focus  of  this  report's  analyses. 
An  extended  treatment  of  each  of  the  parameters  in  Equation  1  is  given  in  Degnan  [1]. 


2.2  JASON  System  Specifications  and  Orbit  Characteristics 

For  this  analysis,  we  assumed  a  typical  TOPEX/POSEIDAN  orbit:  nominally  circular,  with  66.0° 
inclination,  and  an  altitude  over  the  site  of  1,330  km.  The  site  itself  was  assumed  to  be  Capraia 
Island,  latitude:  43°0'  N;  longitude:  9°4.5'  E.  The  expected  performance  of  the  single  retroreflector 
and  the  22  element  array  was  determined  for  Passes  44  and  85  over  this  location. 

For  the  link  analysis,  we  assumed  that  the  transmit/receive  (T/R)  system  has  the  following 
characteristics:  Wavelength:  532  nm  and  1064  nm;  Average  Energy:  100  mJ;  Pulse  width:  100  ps 
at  532  nm;  Transmission  Path  Efficiency:  50%;  Telescope  Diameter:  13  cm;  and  Beam  Divergence 
(Full  angle):  100  ytirad.  We  used  MODTRAN  to  estimate  the  atmospheric  transmission  Ta  for 
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IVODTRAN  AtTiospheric  Absorption 


Figure  2:  MODTRAN  estimates  for  atmospheric  transmission  as  a  function  of  elevation  angle 
6eiev  for  a  clear  sky  at  532  nm  and  1064  nm  above  Capraia  Island. 

Capraia  Island.  The  atmospheric  transmission  is  shown  as  a  function  of  elevation  angle  in  Figure  2 
for  both  532  nm  and  1064  nm.  We  assumed  that  no  cirrus  clouds  were  present  and,  hence,  Tc  =  1. 

The  system  pointing  uncertainty  was  unspecified,  so  we  analyzed  the  link  for  two  different  values 
of  this  parameter.  Based  on  our  experience  on  the  NRL'SSOR  mount  [2]  where  pointing  jitter  was 
on  the  order  of  an  arcsecond  (4.85/irad),  we  estimated  jitter  to  be  on  the  order  of  2  arcsec  for  the 
analysis  in  one  case.  We  also  analyzed  the  array’s  performance  based  on  the  tracking  uncertainty  of 
the  1988  Transportable  Laser  Ranging  System  (TLRS)  upgrade,  which  is  estimated  to  be  5.4  arcsec. 
As  will  be  seen  in  the  results,  the  difference  in  these  two  levels  of  uncertainty  impact  the  predicted 
photoelectron  and  photon  returns. 

The  JASON  specification  stated  that  the  design  shall  provide  a  return  power  of  100  photons 
using  a  Receiver  efficiency  of  20%  in  the  green  and  2-4%  in  the  infrared,  where  these  percentages 
reflect  the  combined  receive  path  and  detector  efficiencies.  If,  in  fact,  that  is  the  case,  then  we 
were  provided  with  the  effective  receiver  efficiencies  to  produce  the  number  of  photoelectrons, 
not  photons.  The  detector  efficiency  is  the  conversion  efficiency  from  photons  to  photoelectrons. 
Therefore,  we  analyzed  the  performance  of  the  array  in  terms  of  photoelectrons.  Ape,  using  20%  for 
532  nm  and  4%  for  1064  nm.  However,  we  also  provide  photons  received  at  the  telescope,  A-^^telei 
for  an  assumed  13  cm  diameter  aperture.  These  numbers  can  be  converted  by  the  product  of  the 
receiver  path  efficiency  and  detector  efficiency,  to  the  number  of  detected  photoelectrons. 

2.3  Target  Diffraction  Effects  and  ctlrcs 

The  typical  target  for  an  SLR  application  uses  one  or  more  retroreflectors.  A  cube  corner  retrore- 
fiector  (or  cube  for  short)  is  made  of  three  mutually  perpendicular  reflective  flats.  The  configuration 
returns  light  along  the  path  of  incidence,  and  therefore,  avoids  a  requirement  for  precise  orientation 
of  the  target  with  respect  to  the  ground  site.  However,  since  the  retroreflector’s  aperture  is  of 
finite  size,  the  reflected  light  is  spread  by  diffraction.  The  transmitted  pulse  travels  hundreds  of 
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kilometers  and  over  the  retrorefiector's  aperture  the  wave  fronts  are  indistinguishable  from  a  plane 
wave.  When  combined  with  the  return  path,  the  geometry  satisfies  the  Fraunhofer  limit.  Conse¬ 
quently,  the  return  pulse  spatial  distribution  is  the  aperture's  Far  Field  Diffraction  Pattern  (FFDP). 
Further,  the  FFDP  follows  standard  diffraction  scaling  in  which  a  smaller  aperture  retrorefiector 
generates  a  broader  FFDP.  The  FFDP's  breadth  is  of  particular  importance  since,  in  combination 
with  the  satellite  velocity  aberration,  it  plays  a  significant  role  in  overall  SLR  system  performance. 

In  Equation  1  the  SLR  target’s  diffraction  characteristics  are  carried  by  the  laser  radar  cross 
section,  (Jlrcs-  Denoting  the  FFDP  by  a{kj:,ky),  the  relationship  between  the  two  quantities  is 


ky) 


47r 


j  Cl  i^kj^ ,  )  I 


(see  Minott,  [3]  [4])  where, 


and 


27T 

kr  —  ~r~  sin  Ojr 

A 


2tt 

ky  =  —SXndy 


(6) 


(7) 


The  angles  9^  and  By  define  the  observation  direction.  If  an  FFDP  is  observed  at  and  t/o  from 
a  range  R,  then  sin  9^  =  Xo/ R  and  sin  9y  =  yo/R  respectively.  One  is  equally  free  to  use  a  position 
and  a  distance,  but  for  general  work  it  is  more  convenient  to  use  angles.  The  results  of  this  report 
are  given  in  terms  of  angular  measure.  Equation  5  assumes  perfect  reflectivity.  A  more  realistic 
description  of  a  real  retrorefiector  is 


(7LRCs{kx,ky)  —  p^\a(kx,ky)f  ,  (8) 

where  p  is  a  reflectivity.  In  the  computations  supporting  this  report,  p  =  0.75  which  is  a  conservative 
estimate  including  manufacturing  errors  as  well  the  reflectivity  loss  itself. 


2.4  Satellite  Velocity  Aberration 

Because  the  satellite  is  moving  with  respect  to  the  ground  station,  the  part  of  the  intensity  distri¬ 
bution  from  the  return  pulse  which  is  observed  at  the  ground  station  is  offset  with  respect  to  the 
central  maximum  of  the  diffraction  pattern.  This  offset  is  called  the  satellite  velocity  aberration  or 
Bradley  effect  [3j.  Figure  3  shows  schematically  the  combined  effects  of  the  diffracted  return  pulse 
and  the  velocity  aberration.  The  velocity  aberration  expressed  as  a  displacement  vector  in  the  far 
field  diffraction  pattern  (or  (Tlrcs)  is  given  by 


k  = 


(9) 


where  vj_  is  the  satellite’s  apparent  perpendicular  velocity.  The  displacement  vector  k  selects  which 
value  of  the  instantaneous  ctlrcs  from  Equation  8  is  used  in  Equation  1  at  each  point  in  time  during 
a  satellite’s  pass  over  a  ground  station. 


Figure  3:  The  SLR  target  size  determines  the  spatial  distribution  of  the  diffracted  return  pulse. 
The  satellite’s  apparent  orbital  velocity  determines  offset  between  the  central  maximum  of  the 
diffraction  pattern  and  the  ground  station. 
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3  SATELLITE  LASER  RANGING  WITH  SINGLE  RETRO  REFLECTOR  TAR¬ 
GETS 

This  section  discusses  the  received  signal  strength  from  a  single  retroreflector  on  a  satellite.  The 
section  first  covers  the  results  from  our  numerical  computation  of  (Tlrcs  fo'"  four  cases  beginning 
with  the  highest  symmetry:  1.  normal  incidence  without  bevels,  2.  tilted  incidence  without  bevels, 

3.  normal  incidence  with  bevels,  and  4.  tilted  incidence  with  bevels.  Each  case  serves  as  a  limit  test 
case  for  the  numerical  routines  of  the  succeeding  cases.  The  second  part  of  the  section  discusses 
the  combination  of  the  ctlrcs  values  with  rest  of  the  SLR  system  to  predict  system  performance. 

3.1  Circular  Retroreflector  Far  Field  Diffraction  Patterns  and  Laser  Radar  Cross 
Sections 

3.1.1  Normal  Incidence 

When  a  plane  wave  pulse  hits  the  retroreflector,  it  is  exactly  reversed  in  direction  along  the  path 
of  incidence.  .4s  discussed  in  Section  2.3,  the  retroreflector’s  finite  aperture  modifies  the  incidence 
plane  wave  pulse  and  returns  a  diffracted  pulse.  At  the  ground  station  the  return  pulse  has  the 
spatial  distribution  of  the  aperture’s  far  field  diffraction  pattern.  A  single  circular  retroreflector 
without  bevel  losses  observed  at  normal  incidence  provides  a  geometry  with  sufficiently  high  sym¬ 
metry  that  an  analytic  expression  exists  for  the  FFDP.  The  FFDP  is  the  .Airy  function  (Born  and 
Wolfe,  [-5];  Goodman,  [6]),  and  hence. 


c^LRCs(^-)  = 


Air. A? 


(10) 


where  A  is  the  area  of  the  retroreflector,  r  is  the  retroreflector  radius,  and 

^  sin  0  .  (11) 

.At  normal  incidence,  <7lrcs  azimuthal  symmetry  and  Equation  10  is  a  function  of  the  single 
magnitude  variable  k  only. 

Figure  4  shows  crrcs  s-s  function  of  k  (a  radial  slice  in  A:^,  ky)  for  2.54  cm  and  1  cm  aperture  sizes. 
While  both  curves  were  computed  numerically  with  the  FFDP  routines  described  in  Section  A.l, 
they  agree  with  Equation  10  to  machine  precision.  The  correct  geometry  required  for  a  normal 
incidence  observation  can  occur  only  when  the  satellite  is  directly  above  the  ground  station  and  the 
retroreflector  normal  is  aligned  with  the  nadir  direction.  The  velocity  aberration  for  a  1,330  km 
circular  orbit  at  zenith  is  approximately  k  of  50  //radians.  Thus,  as  can  be  seen  in  Figure  4,  although 
(Tlrcs  for  the  2..54  cm  retroreflector  is  significantly  larger  at  /c  =  0,  the  broader  FFDP  of  the  1  cm 
retroreflector  has  an  equivalent  (Tlrcs  3-t  ^  of  50  //radians  and  greater  constancy  (less  variation) 
throughout  the  pattern.  Figure  5  shows  a  contour  plot  of  (tlrcs  of  a  1  cm  circular  retroreflector  as 
a  function  of  and  ky.  Note  that  the  velocity  aberration  given  by  Equation  9  selects  a  single  kx, 
ky  point  in  Figure  5  for  which  the  ctlrcs  value  would  be  used  in  Equation  1. 


3.1.2  Off-axis  Incidence 

At  all  other  points  in  an  orbital  pass,  the  retroreflector  tips  away  from  normal  incidence,  so  its 
projected  aperture  changes  in  shape  and  decreases  in  area.  Figure  6  shows  how  both  the  shape  of  the 
projected  aperture  and  area  changes  as  a  function  of  incidence  angle  Bind-  Note  that  the  projected 
aperture  decreases  in  size  along  both  axes.  The  decrease  in  projected  aperture  size  increases  the 
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Effective  Retroreflector  Area 


Figure  6:  Top  Panel:  Retroreflector  Aperture  Shape:  1.  Normal  incidence,  2.  20°  from  normal 
incidence,  3.  40°  from  normal  incidence.  Bottom  Panel:  The  normalized  retroreflector  aperture 
area  as  a  function  of  incidence  angle  Oind' 


angular  extent  of  the  FFDP  proportionally  and  decreases  its  overall  magnitude.  The  change  in 
aperture  shape  eliminates  azimuthal  symmetry  so  an  analytic  closed  form  expression  for  the  FFDP 
(e.g.  Airy  function)  is  no  longer  possible.  Appendix  A.l  describes  a  method  to  numerically  compute 
the  FFDP. 

Figure  7  shows  ctlrcs  of  a  tilted  1  cm  circular  retroreflector  computed  numerically  without 
bevel  losses-  The  azimuthal  symmetry  present  for  normal  incidence  has  been  reduced  to  a  two-fold 
rotational  symmetry  as  seen  from  the  ellipticity  of  the  contours.  The  central  maximum  around 
=  ky  =  0  remains  smoothly  rounded.  However,  there  is  a  sizable  variation  in  aLRcs  one  holds 
the  magnitude  of  k  fixed  and  sweeps  around  in  ky.  Physically  this  corresponds  to  holding  the 
retroreflector  axis  at  a  fixed  tilt  while  sweeping  the  observation  direction  around  at  a  fixed  angle. 
This  variation  is  shown  in  Figure  8  for  a  10°  tilt  from  normal  incidence  and  |k|  of  50  /^radians  as 
a  function  of  observation  direction  angle. 

3.1.3  Bevel  Los.ses 

Beveled  edges  are  required  at  the  reflecting  face  plane  intersections  during  the  polishing  of  real 
retroreflector  cubes.  Adding  bevels  introduces  narrow  loss  regions  within  the  aperture.  More 
general  versions  of  our  numerical  simulations  include  these  losses  in  both  normal  incidence  or  off- 
axis  cases  (Sections  A. 1.3  and  A. 1.4).  At  normal  incidence,  the  three  bevel  loss  regions  and  their 
reflections  generate  a  six-fold  rotation  symmetry  in  diffraction  pattern  and  <7lrcS‘  Since  the  bevels 
are  relatively  narrow,  the  six-fold  symmetry  is  more  pronounced  at  larger  values  of  |kl  beyond 
the  first  diffraction  minimum.  Figure  9  shows  o*lrcs  normal  incidence  of  circular  retroreflector 
computed  numerically  with  losses  from  0.008"  bevels.  The  variation  in  c^lrcs  as  one  sweeps  around 
in  kx,  ky  is  shown  in  Figure  10  for  0.008"  bevels  at  |k|  of  SO/xradians 

As  can  be  seen  from  the  comparison  of  Figures  8  and  10,  the  effect  of  tilting  the  retroreflector  is 
much  more  significant  than  adding  the  bevels.  While  our  numerical  computation  treats  the  bevel  s 
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Figure  7:  Plot  of  (Tlrcs  of  a  1  cm  circular  retroreflector  as  a  function  of  and  ky  given  in  terms 
of  /iradians,  with  10°  tilt.  Ellipticity  of  the  contours  is  the  result. 


1  cm  retroreflectors,  50  microradian,  10  degree  tilt,  no  bevel 


Figure  8:  Variation  in  (Tlrcs  at  10°  off-axis  from  normal  incidence  of  a  1  cm  circular  retroreflector 
at  |k|  of  50  /rradians.  The  variation  in  cjlrcs  from  its  average  caused  by  a  fixed  10°  tilt  by  a 
rotation  of  the  observation  direction  is  42%.  .4t  the  instant  of  pulse  reflection,  asatellite  s  apparent 
perpendicular  velocity  determines  the  observation  direction. 
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Figure  9:  Plot  of  ctlrcs  of  a  1  cm  circular  retrorefiector  with  0.02032  cm  (or  0.008")  bevels  as  a 
function  of  kx  and  ky  given  in  terms  of  /xradians  for  normal  incidence. 


1  cm  retroreflectors,  50  microradian,  0  degree  tilt,  8  mil  bevel 


Figure  10:  The  variation  in  (Tlrcs  at  normal  incidence  of  a  1  cm  circular  retrorefiector  at  jk]  of 
50  /iradians  with  0.02032  cm  (or  0.008")  bevels.  The  variation  from  the  average  ctlrcs  caused  by 
the  bevels  is  0.25%. 
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Figure  11:  ctlrcs  of  a  1  cm  circular  retroreflector  (top  panel)  and  0siev  (bottom  panel)  as  a  function 
of  time  during  Pass  44.  At  PCA,  OeUv  =  89.38°  so  this  is  effectively  a  zenith  pass.  Note  the  impact 
of  retroreflector  Field  of  View  (FOV)  limit  on  ctlrcs- 


effect  exactly,  to  first  order,  the  bevels  could  also  be  considered  as  simply  a  scalar  loss  of  reflecting 
area  combined  with  a  slightly  reduced  field  of  view.  The  combination  of  both  retroreflector  tilt 
and  bevel  losses  requires  a  slightly  more  complex  numerical  simulation,  but  does  not  introduce  any 
qualitatively  new  features  to  the  FFDP  or  (Tlrcs- 

3.2  Single  Retroreflector  Performance  in  Orbit 

The  local  geometry  at  the  satellite  and  apparent  velocity  for  a  given  ground  station  must  be  used 
to  drive  the  (Tlrcs  routine  to  compute  the  performance  expected  for  a  single  retroreflector  on  orbit. 
The  range  vector,  apparent  velocity  vector,  and  local  elevation  angle  from  nadir  are  sufficient  to 
determine  the  kx,  ky  point  and  what  form  of  ctlrcs  wffi  occur  at  a  given  point  in  time  during  a 
satellite’s  orbital  pass.  (See  Section  A.2  for  detailed  development.)  These  two  vectors  and  angle 
were  computed  as  a  function  of  time  for  all  passes  during  an  ten  (9.92)  day  ground  repeat  track  for 
Point  of  Closest  Approach  (PC.A)  Oeiev  >  These  passes  were  provided  by  NRL  Code  8233  for 
the  TOPEX/Poseidon  spacecraft  over  Capraia  Island  is  specified  in  Section  2.2.  The  (Tlrcs  values 
for  a  single  1  cm  nadir-facing  retroreflector  with  0.008"  bevels  are  shown  in  Figure  11  for  Pass  44. 
As  can  be  seen  from  the  figure,  the  single  retroreflector’s  Field  of  View  (FOV)  cutoff  (56.7221°)  is 
a  dominant  constraint;  (tlrcs  is  held  to  zero  during  those  periods  in  the  pass  when  <  38  . 

The  number  of  detected  photons,  AX, tele,  and  photoelectrons,  Ape  for  a  single  1  cm  retroreflector 
with  0.008"  bevels  for  Pass  44  are  shown  in  Figures  12  and  13.  The  performance  is  evaluated  based 
on  both  2  arcsec  and  5.4  arcsec  pointing  uncertainties  and  50%  receiver  path  efficiency  at  532  nm 
as  discussed  in  Section  2.2.  As  can  be  seen  from  the  Figure  13,  detection  of  5  photoelectrons  may 
not  be  possible  below  of  57°.  Thus,  orbital  sampling  will  be  severely  impacted  over  the  site 
for  a  single  1  cm  retroreflector. 
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Time  (sec) 

Figure  12:  A-^,teie  for  a  13  cm  diameter  aperture  and  1  cm  circular  retroreflector  as  a  function  of 
time  during  Pass  44  assuming  either  a  2  arcsec  or  5.4  arcsec  pointing  uncertainty. 
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Time  (sec) 


Figure  13:  A"pe  of  a  1  cm  circular  retroreflector  as  a  function  of  time  during  Pass  44  assuming  either 
a  2  arcsec  or  5.4  arcsec  pointing  uncertainty. 
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Figure  14:  The  NRL  LEO  retroreflector  array  using  22  individual  retroreflectors  on  the  surface  of 
hemisphere.  This  array  has  a  FOV  of  108°  and  (Tlpcs  ^  for  ^  20  . 

4  RETROREFLECTOR  ARRAY  DESIGNED  FOR  LEO 

In  order  to  maximize  opportunities  for  ranging  and  orbital  sampling  over  the  site,  a  retroreflector 
array  can  be  designed  to  provide  a  robust  return  margin  from  horizon-to-horizon,  if  required.  The 
NRL  LEO  retroreflector  array  was,  in  fact,  designed  to  provide  a  ctlrcs  >  10'‘  m^,  through  the  arc 
defined  by  ±20°  for  a  given  pass.  This  array  is  comprised  of  22  1  cm  retroreflectors  specifically 
oriented  in  a  hemisphere.  The  retroreflectors  are  inset  slightly  which  results  in  some  vignetting 
{i.e.  c  /  0  in  Equation  .A.-8).  A  photograph  of  the  array  is  shown  in  Figure  14. 

The  array’s  multiple  retroreflectors  are  specifically  placed  to  provide  a  Field-of-View^  of  108°, 
yet  not  impact  timing  (see  Section  4.3).  A  mechanical  drawing  of  the  array  is  shown  in  Figure  15. 
The  retroreflectors  are  in  three  rings.  The  first  ring  is  at  16°  with  four  retroreflectors  placed  at  0°, 
90°,  180°,  and  270°,  in  azimuth.  The  second  is  at  32°  with  eight  retroreflectors  placed  at  22.5°, 
67.5°,  112.5°,  157.5°,  202.5°,  247.-5°,  292.5°,  and  337.5°  in  azimuth.  Finally,  the  third  is  at  48°  with 
ten  retroreflectors  placed  at  18°,  54°,  90°,  126°,  162°,  198°,  234°,  270°,  306°,  and  342°  in  azimuth. 

.As  with  the  single  retroreflector,  far  field  diffraction  effects  and  velocity  aberration  must  be 
taken  into  account  to  compute  predicted  performance.  However,  as  will  be  shown,  returns  from 
multiple  retroreflectors  for  a  given  pulse  mitigates  nulls  in  the  (Tlrcs  and  permits  a  much  larger 
FOV. 


4.1  ctrrcs  of  the  Retroreflector  Array 

The  return  pulse  is  the  sum  of  the  contributions  from  several  retroreflectors  in  the  array.  The  total 
FFDP  a  for  the  array  of  L  retroreflectors  at  a  given  instant  of  time  is 

L 

fci 


(12) 
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where  the  a/s  are  the  phase  angles  for  each  retroreflector  for  a  given  incidence  direction.  From 
Equation  8,  the  instantaneous  laser  radar  cross  section  is 


Upon  substitution: 


f^LRCS 


(13) 


(14) 


(15) 


During  the  brief  interval  r  while  the  laser  pulse  is  reflecting  from  the  array,  small  variations 
in  orientation  will  cause  essentially  random  fluctuations  among  the  phase  angles  for  the  individual 
retroreflectors.  Consequently,  a  time-averaged  LRCS  will  simplify  to: 


<7LRCS 


(16) 


(IT) 


For  SLR  measurements  with  precisions  significantly  tighter  than  ±1  cm,  r  may  be  become  small 
enough  that  the  phase  integral  can  not  be  considered  as  a  Kronecker  delta  For  JASON, 

however,  using  a  random  phase  approximation  is  sufficient,  and  the  array  s  laser  radar  cross  section 
can  be  evaluated  using  the  incoherent  sum  of  the  contributions  in  Equation  li. 

Equation  17  is  a  convenient  form  for  numerical  evaluation.  The  routines  developed  and  tested 
in  the  single  retroreflector  can  be  used  directly  to  compute  the  fimS.  The  primary  addition  is 
a  routine  which  combines  an  incidence  direction  with  respect  to  the  array’s  axis  with  individual 
retroreflector  directions  to  determine  which  retroreflectors  are  illuminated  within  their  FO\'  limits 
and  at  what  angles.  It  is  then  straightforward  to  compute  the  relevant  a^s  and  to  sum  to  obtain 
ctlrcs-  Figure  16  shows  <tlrcs  of  the  array  (using  retroreflectors  with  0.008"  bevels)  as  computed 
numerically  when  viewed  along  the  array’s  axis. 


4.2  Array  Performance  in  Orbit 

.As  with  the  single  retroreflector,  local  geometry  and  velocity  aberration  combine  to  determine 
o'LRCS  ^t  any  given  instant.  Similarly,  Equation  1  can  be  used  to  determine  Ape  as  a  function  of 
time  during  a  pass  over  a  given  ground  site.  As  in  the  case  of  the  single  retroreflector,  the  orbital 
parameters  for  the  TOPEX/Poseidon  spacecraft  over  an  ten  day  ground  repeat  track  over  Capraia 
Island  were  used  to  characterize  the  array's  performance.  For  an  array,  the  local  angle  from  nadir 
is  the  angle  between  the  incidence  direction  and  the  array’s  axis.  The  apparent  velocity  vector  is 
similarly  common  to  the  entire  array.  .As  noted  in  the  previous  section,  this  information  is  first 
converted  into  the  coordinate  system  of  each  illuminated  retroreflector  in  the  array  and  then  used 
to  compute  ctlrgs- 

Figures  17,  18,  and  19  show  graphs  of  CLFiCS  and  OeUv  versus  time  for  Pass  44,  Pass  85  and  a 
low  elevation  off-zenith  pass  respectively.  As  can  be  seen  from  Figures  li  and  18,  although  there 
is  a  drop  at  nadir,  ctlrcs  remains  greater  than  10’*  m^  throughout  the  entire  pass.  By  comparison, 
the  response  of  a  single  retroreflector  (Figure  11)  shows  vignetting  by  the  FOV  and  severely  limits 
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Figure  16:  NRL  LEO  Retroreflector  array  log((7LRcs)  with  0.008"  bevels  and  on-axis  incidence. 


Pass  44  —  89.38  PCA  September  11,  1996  12:23 


Figure  17:  ulrcs  of  the  22  element  NRL  LEO  retroreflector  array  (top  panel)  and  0elev  (bottom 
panel)  as  a  function  of  time  for  Pass  44,  at  PCA,  6siev  =  89.38°. 
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Time  (sec) 

Figure  18:  ctlrcs  of  th®  22  element  NRL  LEO  retroreflector  array  (top  panel)  and  ^^lev  (bottom 
panel)  as  a  function  of  time  for  Pass  85.  at  PCA,  =  89.31°. 


Pass  ??  —  30.54  PCA  September  18,  1996  09:23 
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Figure  19:  ctlrcs  of  the  22  element  .\RL  LEO  retroreflector  array  (top  panel)  and  ^elev  (bottom 
panel)  as  a  function  of  time  for  a  low  elevation  pass,  at  PCA,  Oghv  =  30.54°. 
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Figure  20;  A%,teie  for  ^  13  cm  diameter  aperture  and  the  NRL  LEO  retroreflector  array  for  Pass  44 
(top  panel)  and  ^elev  (bottom  panel)  as  a  function  of  time  during  a  pass.  A  pointing  precision  of 
2  arcsec  at  532  nm  was  used  in  the  computations. 


viable  returns  outside  a  ±57°  window.  For  the  low  elevation  pass,  Figure  19  shows  ctlrcs  ±  10^ 
occurring  even  at  the  lowest  elevation  angles. 

Performance  on  orbit  in  terms  of  photon  and  photoelectron  returns  are  summarized  for  dif¬ 
ferent  parameter  sets  in  Figures  20  through  23.  In  these  computations,  we  used  the  MODTRAN 
atmospheric  transmission  values  shown  in  Figure  2.  For  Pass  44,  assuming  a  tracking  uncertainty 
of  2  arcsec,  we  see  from  Figures  20  and  21,  that  we  can  detect  18  photons  and  7  photoelectrons  as 
low  as  24.6°  elevation  at  532  nm.  It  is  of  interest  to  note  that  for  a  possibly  more  realistic  pointing 
accuracy  of  5.4  arcsecond,  minimum  elevation  is  28.3°  to  detect  a  comparable  number  of  photons 
and  photoelectrons.  These  results  can  be  seen  in  Figures  22  and  23. 

Performance  for  Pass  85  was  also  evaluated.  Results  were  computed  for  returns  at  532  nm  and 
1064  nm.  Photon  and  photoelectron  returns  through  the  pass  are  shown  in  Figures  24  and  25.  As 
can  be  seen  in  these  graphs,  the  4%  efficiency  of  the  detector  at  1064  nm  significantly  impacts  the 
arc  sampled.  Predicted  returns  at  5.4  arcsec  pointing  uncertainty  are  shown  in  Figures  26  and  27. 


4.3  Optical  Phase  Center  and  Timing  Precision 

The  orbital  estimation  precision  using  the  SLR  data  type  is  on  the  order  of  a  ±1  cm,  RMS. 
Therefore,  knowledge  of  the  offset  of  the  retro  reflector's  location  with  respect  to  the  spacecraft 
Center  of  Mass  (COM)  is  important.  This  fixed  offset  is  called  the  range  correction  and  will  exist 
whether  a  single  retroreflector  or  an  array  is  used. 

In  the  case  of  an  array,  the  timing  precision  of  the  return  signal  can  also  be  affected  if  more 
than  one  retroreflector  is  illuminated  by  a  pulse.  The  placement  of  retroreflectors  defines  an  Optical 
Phase  Center  (OPC)  as  the  common  origin  of  the  all  the  retroreflector  outward  normals.  To  the 
extent  the  array  retroreflectors  lacks  a  common  OPC,  a  direction-dependent  offset  will  be  introduced 
into  the  system.  Consequently,  to  achieve  centimeter-level  precision  and  accuracy,  the  retroreflectors 
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Figure  21:  A’pe  of  the  NRL  LEO  retroreflector  array  for  Pass  44  (top  panel)  and  ^eiev  (bottonn 
panel)  as  a  function  of  time  during  a  pass.  A  pointing  precision  of  2  arcsec  at  532  nm  was  used  in 
the  computations. 
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Figure  22:  AL.tele  for  a  13  cm  diameter  aperture  and  the  NRL  LEO  retroreflector  array  for  Pass  44 
(top  panel)  and  6'eiev  (bottom  panel)  as  a  function  of  time  during  a  pass.  A  pointing  precision  of 
5.4  arcsec  at  532  nm  was  used  in  the  computations. 
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Pass  44  --  89.38  PCA  September  11,  1996  12:23 
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Figure  23:  Ape  of  the  NRL  LEO  retroreflector  array  for  Pass  44  (top  panel)  and  61eiev  (bottom 
panel)  as  a  function  of  time  during  a  pass.  A  pointing  precision  of  5.4  arcsec  at  -532  nm  was  used 
in  the  computations. 


Pass  85  --  89.31  PCA  September  13,  1996  03:27 
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Figure  24:  for  a  13  cm  diameter  aperture  and  the  NRL  LEO  retroreflector  array  for  Pass  85 

(top  panel)  and  ^eiev  (bottom  panel)  as  a  function  of  time  during  a  pass.  A  pointing  precision  of 
2  arcsec  was  used  in  the  computations  at  both  532  nm  and  1.064  pm. 
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Pass  85  --  89.31  PCA  September  13,  1996  03:27 


Figure  25:  Ap^  of  the  NRL  LEO  retroreflector  array  for  Pass  85  (top  panel)  and  (9eiev  (bottom  panel) 
as  a  function  of  time  during  a  pass.  A  pointing  precision  of  2  arcsec  was  used  in  the  computations 
at  both  532  nm  and  1.064  /rm. 


Figure  26:  A^.tele  for  a  13  cm  diameter  aperture  and  the  XRL  LEO  retroreflector  array  for  Pass  85 
(top  panel)  and  ^elev  (bottom  panel)  as  a  function  of  time  during  a  pass.  A  pointing  precision  of 
5.4  arcsec  was  used  in  the  computations  at  both  532  nm  and  1.064  /rm. 


JASOySLR  Target  AnaUisis  vl.O  NRl  Memo  1911 


23 


Pass  85  —  89.31  PCA  September  13,  1996  03:27 


Figure  27:  A'pe  of  the  NRL  LEO  retrorefiector  array  for  Pass  85  (top  panel)  and  ^eiev  (bottom 
panel)  as  a  function  of  time  during  a  pass.  A  pointing  precision  of  5.4  arcsec  was  used  in  the 
computations  at  both  532  nm  and  1.064  jum. 


in  a  multi-element  array  must  be  aligned  to  a  comparable  tolerance  to  a  common  surface. 

The  OPC  of  the  NRL  LEO  array  is  centered  on  the  vertical  axis  6  mm  inside  the  array.  Cen¬ 
tered  about  this  phase  center  is  a  surface  of  reflection.  The  advantage  of  a  hemisphere  of  many 
small  reflectors  is  for  any  angle  there  is  a  retrorefiector  in  close  proximity.  This  condition  allows 
a  relatively  smooth  surface  of  reflection  overall.  The  radius  of  this  surface  was  determined  by  ex¬ 
amining  585  angles  and  calculating  the  point  on  this  surface  for  each  angle.  The  surface  can  be 
approximated  by  a  sphere  with  a  radius  of  2.505  cm.  The  standard  deviation  from  this  sphere  is 
only  0.086  cm.  This  radius  is  the  largest  contributor  to  the  range  correction. 

The  range  correction  also  must  consider  any  pulse  broadening  due  to  the  multiple  retroreflectors. 
Using  small  retroreflectors  allows  for  a  very  compact  arrangement.  The  small  size  limits  the  possible 
pulse  broadening.  The  broadening  results  from  retroreflectors  at  different  locations  on  the  array 
having  a  delayed  return.  For  example,  using  the  nadir  direction  the  ring  of  4  retroreflectors  in  the 
center  are  closer  to  the  ground  then  the  middle  ring  of  8  retroreflectors.  The  time  delay  is  66  ps. 
But  the  return  from  the  middle  ring  is  smaller  due  to  the  larger  angle.  Combing  the  two  returns 
forms  a  broaden  pulse.  The  ring  of  4  has  an  average  ctlrcs  of  o2,  000  m^  and  the  ring  of  8  has 
an  average  ctlrcs  of  only  16,000  m^  for  nadir.  Summing  the  two  returns  for  an  incident  pulse  of 
100  ps  FWHM  is  showm  in  Figure  28.  The  pulse  is  broadened  to  117  ps  FWHM  with  a  7  ps  (1  mm 
)  delay  in  the  peak  amplitude.  The  extra  17  ps  of  width  contributes  to  the  range  error.  A  1  mm 
tw^o  w^ay  range  correction  is  required  to  correct  for  pulse  broadening.  Longer  pulses  w^ould  have 
smaller  effects. 

The  range  correction  error  for  the  time  of  flight  is  the  combination  of  the  size  of  the  sphere 
error  and  the  pulse  broadening.  The  increase  in  the  total  standard  deviation  error  of  the  two  is 
±0.86  mm  and  ±2.6  mm.  Using  Gaussian  statistics  results  in  a  total  error  of  ±2.74  mm.  This  shows 
that  the  array  is  wdth  in  a  ±  centimeter  specification.  A  more  accurate  error  analysis  wdll  require 
detailed  information  including  the  nonlinear  behavior  of  the  receiver  detector  and  electronics. 

It  should  be  noted  that  multiple  retroreflectors  interfere  due  to  the  coherence  of  the  laser  beam. 
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A^ray  Pulse  Timing  Effects 


Figure  28:  Array  Pulse  Timing  Effects 

This  effect  is  termed  satellite  scintillation.  An  incoherent  sum  over  the  array’s  retroreflectors 
generates  the  average  return  signal,  but  a  coherent  sum  is  required  for  determining  the  standard 
deviation  of  the  return  signal.  Having  more  than  two  retroreflectors  in  the  FOV  reduces  the 
deviation.  Since  the  location  of  each  retroreflector  at  a  given  moment  in  time  is  not  known, 
random  phases  are  assigned  to  each  retroreflector.  This  is  equivalent  to  moving  each  retroreflector 
less  than  l/rm.  This  was  done  analytically  for  the  (Tlrcs  each  point  in  the  pass  and  was  repeated 
for  100  random  phases  for  each  retroreflector.  The  resulting  standard  deviation  of  the  LRCS  was 
less  than  12%  of  the  average  return  signal  for  all  855  samples  of  Pass  44.  See  Figure  29.  This 
deviation  is  small  compared  to  atmospheric  scintillation. 

4.4  Further  Comparison  of  Single  Retroreflectors  and  Retroreflector  Arrays 

Using  an  array  eases  the  minimum  pass  geometry  requirements  significantly.  As  has  been  shown, 
an  array  can  provide  viable  returns  at  much  lower  elevation  angles.  Using  the  large  set  of  pass 
geometries  from  the  ten  day  repeat  cycle  as  a  sample  set,  the  number  of  passes  that  exceed  a  given 
value  of  Beiev  at  PCA  as  a  function  Beiev  is  shown  in  Figure  30.  In  other  words.  Figure  30  shows 
the  number  of  passes  available  versus  ^ekv  at  PCA.  From  Figure  11  the  single  retroreflector  (Tlrcs 
places  limit  of  B^iev  >~  40°.  For  an  array,  additional  tracking  passes  can  be  obtained  for  a  given 
length  of  observation  time  and,  consequently  a  better  orbit  solution  becomes  available. 
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5  SPACE  QUALIFICATION  OF  THE  NRL  LEO  RETROREFLECTOR  AR¬ 
RAY 

The  NRL  LEO  Retroreflector  Array  was  qualified  for  space  flight  with  a  number  of  tests.  The  Qual¬ 
ification  Unit  is  identical  to  the  flight  articles  except  that  the  optics  are  slightly  flawed.  Specifically, 
the  bevels  were  below  the  0.004"  specification  on  the  retroreflectors.  The  flight  qualified  retrore- 
flectors  were  qualified  separately  by  Zygo. 

The  data  taken  from  the  tests  described  below  was  analyzed.  No  thermal-related  problems  were 
detected  and  no  negative  margins  of  safety  for  strength  or  fatigue  were  found.  A  more  complete 
description  of  the  qualification  results  can  be  found  in  Appendix  B. 

The  unit  weighs  0.782  ounces,  including  glass,  mount  assembly,  and  bolts  for  mounting. 

5.1  Random  Vibration  Tests 

To  conduct  the  random  vibration  tests,  the  unit  was  mounted  on  an  electrodynamic  shaker.  The 
axial  and  lateral  axes  were  shaken  independently'.  Each  vibration  was  2  minutes  in  duration.  The 
axial  level  was  13.92  G.  rms.  The  lateral  level  was  13.92  G,  rms.  as  well. 


5.2  Thermal  Vacuum  Tests 

The  unit  was  mounted  on  a  thermally'-controlled  heat  sink  and  thermocouples  were  used  to  monitor 
and  control  temperature.  The  vacuum  was  less  than  0.00001  Torr.  The  unit  was  tested  through  6 
cycles  from  —100®  C  to  71°  C. 

5.3  Pyroshock  Tests 

The  unit  was  mounted  on  an  electrodynamic  shaker  and  again  the  two  axes  where  shocked  inde¬ 
pendently  in  axial  and  lateral  dimensions.  The  peak  levels  of  shock  were  6,000  G’s. 
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6  CONCLUSIONS  AND  CONCERNS 

Based  on  the  analysis  presented  in  this  report,  we  predict  that  JPL  and  JASON  will  be  able  to 
conduct  laser  ranging  successfully  using  the  NRL  LEO  array  and  the  FTLRS  system  at  532  nm,  day 
and  night  at  the  TOPEX  orbit.  However,  closing  the  link  in  the  infrared  using  the  specifications 
provided,  particularly  in  the  daytime,  will  require  photon  counting  for  the  small  FTLRS  aperture. 
Furthermore,  use  of  this  telescope  may  preclude  laser  ranging  at  the  lower  elevation  angles. 
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A  DIFFRACTION  PATTERN  AND  PASS  GEOMETRY  CALCULATIONS 

A.l  Circular  Retro  reflector  Far  Field  Diffraction  Patterns 

A  convenient  method  to  develop  the  analytic  basis  for  the  numerical  determination  of  far  field 
diffraction  pattern  begins  with  the  simplest  case  and  then  adds  refinements.  With  this  approach, 
numerical  implementation  of  each  stage  can  serve  as  a  limiting  test  case  for  the  next  stage.  Our 
starting  case  is  normal  incidence  wdthout  bevel  losses  and  we  will  finish  wdth  tilted  incidence  wdth 
bevel  losses.  As  has  been  noted  in  the  body  of  this  report,  introducing  a  tilt  in  the  incidence 
direction  is  much  larger  effect  than  are  the  bevel  losses. 


A.  1.1  Normal  Incidence  Without  Bevel  Losses 

Consider  a  circular  retroreflector  wdth  bevel  loss  regions  directly  overhead.  The  two-dimensional 
Fourier  transform  of  the  complex  reflectance  of  the  aperture,  effectively  the  retroreflector  s  far  field 
diffraction  pattern,  A{x,y)  is  given  by 


ky) 


=  f  I  dxdy.A{x,y)e  . 

J  J  aperture 


(A-l) 


Using  the  function  c[x)  =  v/r^  -  wTich  defines  the  top  half  of  the  circular  aperture,  the  double 
integral  can  be  rewritten  as 

d{kr.ky)  =  r  r'""' 
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(A-3) 

(A-4) 

(A-5) 

(A-6) 

(A-7) 


When  Equation  A-7  is  integrated  numerically,  the  absolute  magnitude  of  d{kx,ky)  determined,  and 
then  converted  to  a  cross  section,  the  result  is  the  azimuthally  symmetric  Airy  function  of  Figure  4. 
The  numerical  routine  that  implements  this  geometry  serves  as  a  limit  test  case  for  the  tilted  and 
normal  incidence  with  bevel  losses  routines. 


A.  1.2  Tilted  Incidence  Without  Bevel  Losses 

When  a  circular  retroreflector  is  tilted,  the  aperture  changes  shape  to  the  overlap  region  of  two 
offset  circles.  The  tilt  away  from  normal  incidence  can  be  related  to  circles  offset  xq  from  the 
center  of  the  retroreflector.  For  a  retroreflector  of  depth  d  wdth  a  face  recessed  by  distance  c  and 
index  of  refraction  n.,  the  expression  relating  xq  to  the  angle  Oind  betw-een  the  surface  normal  and 
the  incidence  direction  is 


.To  =  csin  Oinci  +  dsin  On  , 


(A-8) 
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where 


Oj.  =  arcsin 


(A-O) 


For  a  circular  retroreflector  of  radius  r,  d  =  \/2r.  The  retroreflector’s  Field  Of  View  (FOV)  is 
limited  by  the  requirement  that  Xq  <  r.  For  a  1  cm  diameter  retroreflector  with  n  =  1.46  and 
c  =  0.1016  cm  (or  0.040"),  the  FOV  cut  off  is  at  Ojnei  =  58.579°.  The  aperture  boundary  for  the 
tilted  retroreflector  is  defined  by  c(x)  =  ±\/ r-  —  (,r  —  .ro)-  for  j  <  0  and  c(x)  = 
for  a;  >  0.  The  aperture’s  shape  is  that  of  a  cat's  iris  or  an  .Vnierican  football  profile  shown  in  the 
upper  panel  of  Figure  6  (page  9). 

Integrating  each  side  of  the  aperture  separately  yields  a  tw’o  dimensional  Fourier  transform  for 
tilted  incidence  with  four  integrals  given  by 


ky^ 
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kti 


f  coskxXs\x\  ikyJr^  -  (.t  -  .r’o)^  )  dx 
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'2i 
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j  sin  k:rxsm  yky\Jr'^  —  (a  -  dx 

h  j  sin  kxX  sin  ^ky^Jr'^  —  [x  +  dx 


(.4-10) 


Equation  A-10  yields  reasonable  looking  transform  and  limits  to  an  Airy  function  as  Xq  -4  0. 
Figure  8  (page  10)  shows  the  sizable  variation  in  the  laser  radar  cross  section  of  a  tilted  circular 
retroreflector  as  one  sweeps  and  ky  with  a  constant  magnitude  of  k  of  50  ^radian  (be.  roughly 
a  satellite’s  velocity  aberration  for  a  1,100  km  high  orbit). 


A.  1.3  Normal  Incidence  With  Bevel  Losses 


Consider  a  circular  retroreflector  with  bevel  loss  regions  directly  overhead.  There  are  six  regions, 
symmetric  about  the  a;-axis,  that  can  be  used  to  compute  the  two-dimensional  Fourier  transform 
for  the  aperture.  These  regions  are  shown  in  Figure  31. 

The  two-dimensional  Fourier  transform  of  the  complex  reflectance  of  the  aperture  A{x,y)  is 
now  given  by 


d[kj..  ky]  — 


(.4-11) 


where 


(.4-12) 
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Figure  31:  On-axis  retroreflector  with  bevel  geometry.  The  six  regions  are  numbered  left  to  right. 
Regions  3  and  4  are  each  made  of  two  pieces. 


and 

l~  i^x)  = -^{x  —  w)  .  (A-13) 

The  two  integration  limits  along  x  determined  by  the  intersection  of  /+  (x)  and  l~[x)  and  the  circular 
aperture  are  given  by  xi  =  |  -H  Vl'2r'^  -  3u’^^  and  X2  =  \  +  '/l2r^  —  3u-^^  respectively. 

Rather  than  working  with  all  six  terms  in  Equation  A-11  at  once,  I'll  reduce  them  to  forms 
convenient  for  nunierical  integration  individually.  The  six  terms  are: 
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Both  ai(kj:,  ky)  and  aaikx,  ky)  have  upper  and  lower  limits  of  the  defined  by  the  circular  aperture 
as  in  Equation  A-2  for  the  integration  along  y.  Hence,  these  two  terms  can  be  immediately  simplified 
to  the  same  form  as  Equation  A-7  yielding 


di{kx,ky)  =  cos  kx'X  sin  (ky  \J  r-  -  dx 

ky  J  —  r 

—  —  /  sin  kj^x  sin  r-  —  dx 
ky  J—r 


(A-20) 


and 
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(A-21) 


Note,  of  course,  that  the  x  integration  limits  in  Equations  A-20  and  A-21  are  different  than  those 
in  Equation  A-7.  Similarly,  d2{kx,ky)  and  £i5(Av,A!,)  can  be  reduced  to  integrals  over  i'  of  the  form 


^2  (^’xi 


2  /I  \ 

=  “7“  /  cos  kjrX  sin  ky—^{x  +  w)  1  dx 

ky  J  —  X2  \  \/3  / 

'2i  /I  \ 

+  —  y  s\n  kj^xsm  yky-^{x  +  w)  j  dx 


(A-22) 


and 
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Note  the  sign  change  of  the  real  and  imaginary  parts  of  Equation  .A-22. 

Simplifying  dsik^^ky)  and  ai{kx,ky)  can  be  done  with  the  same  approach.  For  a^{kx,ky] 
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Finally,  repeating  these  steps  for  as{k:^,ky)  yields  the  similar  form 

^3  ^  ^y) 


r-iul2 

/  e-‘^^Uy  + 

/  e-^^y^dy 

J-Xl  [' 

J  —  c[x) 

J-l-(x) 

~  I  y*  <*08  krX  sin  (k.y  \/ dx 

r-w/2  /  f  \ 

+  j  cosA:j..xsin  -  w)  j  dx 

-iiL 


2i  f  /-“'/2  . 

SI 

r-ujj2 


sin  kxXsin  (^kyV -  x^'j  dx 


+ 


/—IL\ 
-J7l 


sin  kxXsin  (  ky-j={x  —  wO  )  dx 
v3 


:A-29) 


(A-30) 


(A-31) 


The  integrals  in  Equations  A-20,  A-21,  A-22,  A-23,  A-29,  and  A-31  are  to  be  evaluated  nu¬ 
merically.  Figure  10  shows  the  variation  in  the  <tlrcs  of  a  beveled  circular  retroreflector  as  one 
sweeps  kx  and  ky  with  a  constant  magnitude  of  k  of  50  /iradian  (i.e.  roughly  a  satellite  s  velocity 
aberration).  In  comparison  to  the  variation  seen  in  Figure  8  (page  10),  the  bevels  do  not  cause 
that  much  change  in  the  magnitude,  but  do  introduce  a  six-fold  symmetry. 


A.  1.4  Tilted  Incidence  With  Bevel  Losses 

For  a  circular  retroreflector  with  bevel  loss  regions,  off-normal  incidence  angles  tilted  about  the 
y-axis  result  in  an  aperture  that  is  a  combination  of  the  cat  s  eye  boundary  and  bevel  loss  regions. 
The  position  of  the  bevel  loss  regions  is  unchanged  and  Equation  A-8  gives  the  correct  offset  in 
this  case  as  well.  Note  that  the  FOV  limit  now  becomes  Xq  <  r  —  j.  For  a  1  cm  diameter 
retroreflector  with  n  =  1.46,  c  =  0.1016  cm  (or  0.040"),  and  w  =  0.01016  cm  (or  0.004")  the 
FOV  cut  off  is  at  e^nd  =  57.63865°.  Since  the  normal  incidence  case  with  bevel  losses  required 
distinct  integration  regions  on  each  side  of  the  aperture,  the  same  six  regions  used  for  Equation  A- 
11  can  also  be  used  in  this  case.  The  two  integration  limits  along  x  determined  by  the  intersection 
of  /■*■(.?)  and  l~(x)  and  the  aperture  change  to  aq  =  1 1-  ~  5(3'’  ~  *o)^|  and 

X2  =  I  -  ^’o)  +  -  i(u)  +  a-'o)^}  respectively.  The  two  dimensional  Fourier  transform 

for  tilted  incidence  can  be  expressed  in  an  analogous  form  to  Equation  A-11,  where  the  resulting 
six  terms  are  given  by: 

eii[kx^ky']  —  T”  f  cos  sin  iky^Ji  ~  (a  :^o)'  j  ^x 

ky  J—(r  —  Xo)  \  / 

—  —  [  sinA:.rasin  ikyJr'^  -  (a  -  ao)^  )  dx 
ky  y_(r-xo)  V  / 


(A-32) 
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and 


a2(k,r,k,f)  =  I  coskrxsln  (^ky^{x +  dx 

4-  ^  y  ^  sin  A;,,. a: sin  ^ky-^{x  +  w')^  dx  . 
d3(kx,ky)  =  f{/  cosA:^.Tsin  (k.y^Jr^  -  (x  -  .ro)'^^  dx 


+  j  COS sin  (^ky-^{x  -  fia’| 

-  ^  I y  sin  fcj-a;  sin  (^ky\J -  (a'  -  a;o)^^  dx 
+  j  sin  Ava;sin  -  u-)^ 

d4{kx,ky)  =  Y  |y  ^^cosA^a-sin  {ky^r^  -  (a;  +  a;o)'^^  dx 

-  y  cosfcra-sin  dx^ 

—  Y  |y  sin  A’,ra;sin  (^ky^Jr'^  -  (a;  +  a;o)^^  dx 
—  j  sin  kj-x  sin  t/a’| 

do{kx- ky)  =  Y  Ju 

2i  .  /  1  \ 

~  Y  J  sin  A’j.-'i- sin  —  t(-)J  da;  , 


2  /*(^-^o)  /  j - \ 

ae[kj:.ky)  =  —J  cos  kj^x  sin  i^kyyjr^  -  {x  +  xo)‘^  j  dx 


cos  Kj^X  I 

2t  /■('— ^o) 

y  dx2 


—  -r-  f  sin  kj:X  sin  (  ky\lr"^  -  (a;  +  a;o)^  )  dx  . 

ku  Jxo 


Equations  A-32  to  A-37  limit  properly  to  the  untilted  case  as  a;o  ->  0.  (1  think.) 


;A-33) 


(A-34) 


(A-35) 


(A-36) 


(A-37 


A. 2  Converting  Pass  Geometry  to  k^ky  Space 

At  a  given  instant  the  orientation  and  velocity  of  the  satellite  with  respect  to  the  ground  station 
determines  effective  aperture  of  the  retro  reflector  and  the  velocity  aberration  respectively. 

For  simplicity,  consider  a  nadir  pointing  retrorefiector  and  use  a  fixed-earth  coordinate  system. 
Hence,  the  ground  station's  location  vector  is  a  constant  denoted  by  G  with  fixed  spherical  coor¬ 
dinates  of  G,  6g.,,  and  (pg,.  Note  that  9g^  is  the  ground  station's  colatitude  measured  down  from 
the  north  pole  rather  that  the  latitude  from  the  equator.  The  satellite  s  position  is  denoted  by  S 
with  spherical  coordinates  of  S,  and  Og/c-  The  satellite  velocity  v  includes  the  effect  of  ground 
station  motion.  Finally,  the  range  vector  between  satellite  and  ground  station  is  R  =  S  -  G. 
Figure  .32  illustrates  the  geometry  and  the  local  coordinate  system  at  the  satellite. 
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Figure  32:  The  left  panel  shows  the  global  definitions  of  G,  S,  and  R.  The  right  panel  shows  the 
coordinate  system  local  to  the  satellite. 

The  angle  of  incidence  9inci,  which  defines  the  aperture  through  Equation  A-8,  is  the  angle 
between  S,  and  R.  Since  G,  S,  and  R  form  a  triangle,  O^nd  can  be  related  to  the  angular  coordinate_s 
of  G  and  S  and  the  elevation  angle  6dev  of  the  satellite  at  the  ground  site.  The  angle  between  G 
and  S  is  given  by- 

cos  7  =  cos  cos  6gs  +  sin  sin  6gs  co&[(psi^  -  <Pgs)  ■  (A-38) 

Hence, 

^elev  +  ^  +  7  +  ^inci  —  ( A-39) 

or  ^ 

^elev  +  1'  +  ^inci  —  '  (.4-40) 

In  practice,  the  values  of  Bind  and  as  a  function  of  time  during  the  pass  are  computed  directly 
by  orbital  dynamics  software.  When  working  with  a  retroreflector  array,  Bind  is  the  angle  between 
R  and  the  array’s  axis.  Hence,  one  must  do  a  little  further  computation  to  determine  the  angle  of 
incident  for  each  retroreflector  with  respect  to  R. 

The  apparent  velocity  of  the  retroreflector  can  be  expressed  in  terms  of  the  velocity,  the  ground 
station  position  and  satellite  position.  Once  the  velocity  components  in  a  satellite  local  coordinate 
sy'stem  are  determined,  the  values  of  k^.  and  ky  that  defined  the  ground  stations  position  in  the  far 
field  diflfraction  pattern  can  be  obtained  immediately.  The  local  coordinate  system  with  x  in  the 
plane  defined  by  S,  and  R 

f  =  -S  , 

R  X  S 

y  = 

IRX  S 

and 

X  =  y  X  i  .  (.4-43) 


(.4-41) 

(A-42) 
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Hence 


R  X  S 

•  V  . 

lRxS| 

(S  -  G)  X  S  ^ 
|(S  -  G)  X  S| 

S  X  G  _ 

|S  X  G| 


(A-44) 


(A-46) 


Similarly 


Finally, 


and 


X  S 

|Rx  S| 

(A-47) 

’(S  -  G)  X  S 

-  y  _ S  •  V 

( A-4R) 

l(S-G)  X  S| 

A  O  V  , 

S  X  S  X  G  _ 

|S  X  G| 

(A-49) 

/27r\  /2vj-\ 

[tJ  [~J 

(A-50) 

/27r\  f'2vy\ 

(A-51) 

The  numerical  routine  to  determine  kj,  and  ky  from  the  results  provided  by  the  orbital  dynamics 
group  is  a  distinct  module  to  simplify  testing. 
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Figure  33:  Lateral  Accelerometer  Data 


B  CODE  8200  SPACE  QUALIFICATION  REPORT 

B.l  Test  Report  For  The  Retro-Reflector  Array 

August  1996 

Robert  B.  Patterson 

U.S.  Naval  Research  Laboratory 

B.1.1  Random  Vibration  Tests 

The  test  article  (RR-ST-1100-01  assembly)  was  subjected  to  two  random  vibration  tests.  One 
in  the  axial  (Z-axis)  direction  and  one  in  the  lateral  (X-axis)  direction.  Since  the  assembly  is  a 
symmetrical  round  dome,  a  test  in  the  Y-axis  direction  was  considered  unnecessary. 

The  test  article  was  assembled  using  flight  quality  parts  with  the  exception  of  the  optics,  which 
are  identical  to  the  flight  quality  units  in  all  physical  characteristics,  but  have  slight  flaws  in  the 
optical  properties.  The  qualification  dome  (RR-ST-1030-1)  was  not  irradiated,  this  will  have  no 
effect  on  the  fit  of  the  parts  to  one  another  or  the  strength  or  stiffness  of  that  part  or  the  assembly. 

The  lateral  test  was  run  first.  The  accelerometer  data,  from  an  accelerometer  mounted  on  the 
shake  table  close  to  the  test  article,  is  shown  in  Figure  33.  The  test  article  survived  the  random 
vibration  test  in  the  lateral  axis. 

The  axial  test  was  run  next.  The  accelerometer  data  is  shown  in  Figure  34.  The  test  article 
survived  the  random  vibration  test  in  the  axial  axis. 


B.l. 2  Shock  Testing 

The  shock  tests  were  next.  Again  only  an  axial  and  lateral  test  were  needed.  The  test  was  performed 
on  the  same  shake  table  that  was  used  to  perform  the  random  vibration  tests.  The  accelerometer 
data  was  obtained  the  same  way,  with  the  accelerometer  mounted  on  the  test  fixture  close  to  the 
test  article. 
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Figure  34:  Axial  Accelerometer  Data 


The  axial  shock  test  was  completed  first.  The  test  article  passed.  The  accelerometer  data  is 
shown  in  Figure  35. 

The  lateral  shock  test  was  then  run.  The  test  article  passed.  The  accelerometer  data  is  shown 
in  Figure  36 


B.1.3  Thermal  Vacuum  Testing 

The  qualification  assembly  was  placed  in  a  small  vacuum  chamber  and  cycled  between  the  survival 
temperature  extremes  to  verify  that  the  assembly  could  survive  the  environments  it  will  see  while 
attached  to  the  host.  This  was  done  by  bolting  the  assembly  to  a  small  "cold  plate  that  had 
its  temperature  controlled  by  allowing  liquid  nitrogen  to  flow  through  copper  tubing  that  was 
soldered  to  the  backside  of  the  copper  “cold"  plate  or  by  applying  power  to  electrical  heaters  that 
are  attached  to  the  front  side  of  the  plate,  or  some  combination  of  the  two.  The  plate  had  two 
thermocouples  mounted  to  it.  One  was  for  feedback  to  the  controller  used  to  control  the  flow  of 
liquid  nitrogen  and  the  other  was  used  as  the  data  for  the  test. 

The  test  article  was  placed  in  the  vacuum  chamber,  which  was  then  pumped  down  to  a  vacuum 
of  at  least  0.0001  torr.  The  assembly  was  then  subjected  to  temperature  extremes  of  71°  C  for 
the  high  temperature  and  —100°  C  for  the  low  temperature.  The  assembly  was  held  at  each 
temperature  extreme  for  a  one  hour  dwell  period  to  insure  that  the  test  article  stabilized  at  the 
test  temperature.  The  test  consisted  of  6  cold  dwells  alternated  with  6  hot  dwells. 

The  qualification  assembly  passed  the  test  without  suffering  any  failures.  The  test  went  accord¬ 
ing  to  plan,  except  when  the  control  valve  stuck  open  and  the  temperature  went  down  to  —157°  C 
during  the  second  cold  dwell.  Since  the  temperature  stayed  below  the  maximum  of  -100°  C,  this 
did  not  negate  this  part  of  the  test.  The  faulty  valve  was  changed  out  during  the  next  transition 
from  cold  to  hot  dwell  and  the  test  continued.  A  time  versus  temperature  graph  of  the  test  profile 
is  shown  in  Figure  37. 
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Figure  3( :  Thermal  Vacuum  Test  Temperature  Profile 

B.2  Test  Procedure  For  The  Retro-Reflector  Array 

June  19.  1996 

Robert  B.  Patterson 

U.S.  Naval  Research  Laboratory 

B.2.1  Test  Definition 

B.2. 1.1  Scope  This  document  contains  the  requirements  and  procedures  for  the  random  vi¬ 
bration,  shock  and  thermal  testing  to  be  performed  on  the  retroreflector  array  in  order  to  qualify 
it  for  space  flight  on  its  host. 

B.2. 1.2  Test  Objective  The  objective  of  this  testing  is  to  insure  that  the  retroreflector  ar¬ 
ray  is  structurally  sound  and  can  withstand  the  launch  and  post  launch  environments  without 
experiencing  any  structural  failure. 

B.2. 1.3  Pass/Fail  Criteria 

1.  No  parts  (spring,  retaining  ring,  retro  seat  or  retro  optic)  ejected  from  the  assembly. 

2.  No  noticeable  permanent  deformation  around  the  seat  area  at  the  optic  opening  in  the  retrore¬ 
flector  housing. 

3.  No  noticeable  permanent  deformation  around  the  retaining  ring  groove  area  of  the  retrore¬ 
flector  housing. 

B.2. 1.4  Test  Description  Since  the  retro-reflector  assembly  is  a  round  part  that  is  symmet¬ 
rical  about  the  axis  of  rotation,  we  will  only  have  to  test  it  in  one  longitudinal  direction  and  the 
axial  direction.  We  will  mount  the  retro-reflector  assembly  to  the  shaker  table  so  that  it  is  oriented 
to  be  tested  in  the  longitudinal  direction.  We  will  then  run  the  test  article  through  the  13.92  Grms 
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Random  Mbration  test  environment.  We  will  then  "hit"  the  test  article  with  an  extremely  sharp 
impulse  in  order  to  simulate  the  shock  environment  (6000  g)  specified  by  the  host  vehicle.  We  will 
now  remount  the  retro-reflector  assembly  so  that  it  is  oriented  in  the  axial  direction  and  repeat  the 
random  vibration  and  shock  tests. 

The  thermal  test  consists  of  bolting  the  retro-reflector  assembly  to  a  "cold"  plate  (that  also 
has  heaters  installed  on  it),  placing  this  assembly  into  a  vacuum  chamber  and  then  cycling  the 
assembly  through  its  survival  temperature  extremes  six  times. 

B.2.1.5  Test  Location  The  test  will  be  conducted  in  the  Naval  Center  for  Space  Technology’s 
Payload  Checkout  Facility,  Building  A-59  at  the  Naval  Research  Laboratory  in  Washington.  DC. 

B.2.1.6  Responsibilities  The  Spacecraft  Engineering  Department  is  responsible  for  the  overall 
performance  of  the  test,  the  data  reduction,  and  production  of  the  final  test  report. 

B.2.1.7  Test  Log  The  test  conductor  shall  maintain  a  daily  test  log  which  shall  include  a 
detailed  record  of  the  test  progress,  significant  events  and  any  deviations  from  the  test  procedures 
outlined  in  this  document. 

B.2.1.8  Photographic  Coverage  Photographs  will  be  taken  of  the  test  specimens,  test  fix¬ 
tures  and  peripheral  test  equipment.  Each  transducer  location  will  be  marked  and  photographed. 

B.2.2  Test  Conditions 

B.2.2.1  General  The  test  will  be  conducted  in  the  ambient  temperature  and  pressure  in  the 
Payload  Checkout  Facility,  Building  A-59,  Naval  Research  Laboratory. 

B.2.2. 2  Applied  Loads  The  test  articles  shall  be  tested  to  the  levels  described  by  the  host 
vehicle  as  listed  in  NCST  Qual.  Procedures,  .4ppendix  A.  The  qualification  unit  shall  be  tested  to 
qualification  levels  and  the  flight  unit  shall  be  tested  to  acceptance  levels. 

B.2.2.3  Test  Report  A  final  test  report  will  be  prepared  to  document  the  test  results.  This 
report  shall  include  overall  test  objectives,  test  article  and  test  setup  descriptions,  reduced  test 
data,  test  anomalies  and/or  significant  observations  and  a  summary  of  the  test  results. 

B.3  Final  Acceptance  Test  Report  Of  The  Retro- Reflector  Array 

November  1996 

Robert  B.  Patterson 

U.S.  Naval  Research  Laboratory 

B. 3.1  Random  Vibration  Tests 

The  flight  units  (RR-ST- 1100-01  assembly)  were  subjected  to  two  random  vibration  tests  each. 
One  in  the  radial  (Z-axis)  direction  and  one  in  the  lateral  (X-axis)  direction.  Since  the  assembly  is 
a  symmetrical  round  dome,  a  test  in  the  \-axis  direction  was  considered  unnecessary. 

The  flight  units  were  assembled  using  flight  quality  parts  and  were  labeled  assembly  A,  B,  and 

C.  The  flight  units  were  tested  to  acceptance  (flight)  levels  6  dB  below  qualification  for  a  total  of 
6.92  GRMS. 
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Figure  38:  Lateral  Accelerometer  Data,  Acceptance  Level  S/n  -A 

The  lateral  test  was  run  first.  The  accelerometer  data,  from  an  accelerometer  mounted  on  the 
shake  table  close  to  the  test  article,  is  shown  in  Figures  38,  39,  and  40.  The  flight  units  survived 
the  random  vibration  test  in  the  lateral  axis. 

The  axial  tests  were  run  next.  The  accelerometer  data  is  shown  in  Figures  41,  42,  and  43.  The 
flight  units  survived  the  random  vibration  test  in  the  axial  axis. 

B.3.2  Shock  Testing 

The  shock  tests  were  next.  Again  only  an  axial  and  lateral  test  were  needed.  The  test  was  performed 
on  the  same  shake  table  that  was  used  to  perform  the  random  vibration  tests.  The  accelerometer 
data  was  obtained  the  same  way,  with  the  accelerometer  mounted  on  the  test  fixture  close  to  the 
test,  article. 

The  axial  shock  test  was  completed  first.  The  flight  units  passed.  The  accelerometer  data  is 
shown  in  Figures  44,  45,  and  46. 

The  lateral  shock  test  was  then  run.  The  flight  units  passed.  The  accelerometer  data  is  Fig¬ 
ures  47,  48,  and  49. 


JASO^■SLR  Target  Analiisis  vl.O  ^■RL  Mtmo  7911 


"Cit-i-.el:  acjjas  i:aR,i  r:: 

n  nn  l  iliin le5:  /icn-  j.'.-.  Ir 


Til-'  i.niic*  -..1 


.!•  Vi:.? 


Figure  39:  Lateral  Accelerometer  Data,  Acceptance  Level  S/n  -B 


Figure  40:  Lateral  Accelerometer  Data,  Acceptance  Level  S/n  -C 
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Figure  48:  Lateral  Shock  Accelerometer  Data,  Acceptance  Level  S/n  -  B 
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Figure  49:  Lateral  Shock  Accelerometer  Data,  Acceptance  Level  S/n  -  C 

C  ACRONYM  LIST 

APD  Avalanche  PhotoDiode 

COM  Center  of  Mass 

FFDP  Far  Field  Diffraction  F^attern 

FOD  Field  of  Dreams 

FOV  Field  of  View 

FTLRS  Field  Transportable  Laser  Radar  Station 
FWHM  Full  Width  at  Half  Maximum 
JPL  Jet  Propulsion  Laboratory 
LEO  Low  Earth  Orbit 

LRCS  Laser  Radar  Cross-Section  (or  (tlrsc) 

MOBLAS  MOBile  LAser  ranging  System 

OPC  Optical  Phase  Center 

PCA  Point  of  Closest  Approach 

SLR  Satellite  Laser  Ranging 

SOR  Starfire  Optical  Range 

TLRS  Transportable  Laser  Ranging  System 


